Short-circuit boundary conditions in ferroelectric PbTiOa thin films 
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We examine the application of short-circuit electrical boundary conditions in density functional 
theory calculations of ferroelectric thin films. Modeling PbTi03 films with Pt electrodes, we demon- 
strate that under periodic boundary conditions of supercells, short-circuit conditions for the elec- 
trodes are equivalently satisfied in two repeated slab geometries: a PbTiOs /metal superlattice 
geometry and a periodic metal/PbTi03 /metal/vacuum geometry, where the metal is Pt or SrRuOs. 
We discuss the benefits of each geometry in the study of ferroelectricity in thin films. 
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The application of relevant electrostatic boundary con- 
ditions in first-principles simulations of ferroelectric thin 
films is extremely important i^^^ In many studies, the 
calculations are carried out on systems composed of a fer- 
roelectric thin film sandwiched between two metal elec- 
trodes in short-circuit conditions, 6 '- i.e., the electrostatic 
potential through both electrodes is the same. The short- 
circuit boundary conditions are also relevant to most ex- 
periments. 8 

Under supercell periodic boundary conditions in the 
first-principles studies, one way to satisfy the short- 
circuit conditions is to consider a slab geometry consist- 
ing of repeated ferroelectric/metal units, in which the 
metal layers are thick enough so that copies of the fer- 
roelectric film do not interact with each other^ In this 
geometry, illustrated schematically in Fig. the short- 
circuit boundary conditions are directly enforced, with 
contact between the metal electrodes at the top and bot- 
tom of the ferroelectric film through a bulk-like region of 
the metal. In this case, the entire ferroelectric/electrode 
structure can be relaxed by considering the unit cell 
strain as a variational parameter. 

It is also possible to model ferroelectric thin films us- 
ing an isolated capacitor geometry that contains elec- 
trode/ferroelectric/electrode slabs that are separated by 
vacuum, such as in Fig. QJd. This geometry is computa- 
tionally less efficient than the continuous geometry (re- 
quiring more atoms and larger supercells), but it is ben- 
eficial in that there is no artificial constraint on the unit 
cell strains along the surface normal direction, and there- 
fore no need to rely on the strains to find the relaxed 
geometry. By separating the effects of the positive and 
negative ferroelectric surface charges on the metal be- 
havior, an isolated geometry also facilitates the study of 
work functions at the ferroelectric surface^ and the com- 
parison of complete and partial screening in thin metal 
films. Furthermore, this arrangement is useful for study- 
ing the effects of having different top and bottom elec- 
trodes, since it eliminates the influence of an additional 
inter-electrode interface. However, unlike in the contin- 
uous geometry, it is not obvious whether the necessary 
electrostatic boundary conditions can be achieved in an 
isolated structure. 

In this paper, we demonstrate using density-functional 
theory (DFT) calculations that short-circuit boundary 



conditions are satisfied in both continuous and iso- 
lated geometries. We determine the minimum electrode 
thickness in each geometry for ultrathin ferroelectric 
PbTi0 3 (001) films with platinum (Pt) and SrRu0 3 elec- 
trodes, and we demonstrate that, given sufficiently thick 
electrodes, the desired electrostatic boundary conditions 
are satisfied in both geometries. Furthermore, we show 
that with electrodes at or above the minimum thickness, 
the ferroelectric structure of the perovskite films, as well 
as the electronic behavior of both the metallic and per- 
ovskite regions, is indistinguishable in the two geome- 
tries. 

Our calculations were performed using DFT with the 
generalized gradient approximation as implemented in 
the ab initio code dacapo^, with a plane wave cutoff of 
30 Ry and a 4 x 4 x I Monkhorst-Pack fc-point mesb^, and 
all calculations correspond to T—0 K. The PbTiC>3 struc- 
tures considered consist of five atomic layers, three layers 
of TiC>2 separated by PbO planes stacked in the (001) di- 
rection, as illustrated by the schematic in Fig.^ Exper- 
imentally, cubic SrTiC>3 and tetragonal PbTiC>3 have the 
same in-plane lattice constant (3.905 A), so the interface 
is stress free. To model stress-free growth of tetrago- 
nal PbTiC>3 on SrTi03, the in-plane lattice constant was 
fixed to the theoretical equilibrium bulk value of 3.86 A. 
All atomic coordinates were fully relaxed until the forces 
on each atom are less than 0.01 eV/ A. Polarization is 
perpendicular to the interface, which is energetically fa- 
vorable for this choice of in- plane lattice constant i& The 
interfacial Pt atoms are in the energetically-favored po- 
sitions above the TiC>2 oxygens, and additional Pt layers 
continue the (001) lattice. In the continuous geometry, 
periodic copies of the PbTiC>3 film are separated by a 
Pt electrode of Ipt = 3, 5, 7, or 9 atomic layers^, while 
in the isolated geometry, the PbTiC>3 film is sandwiched 
between two Pt electrodes, each of thickness /pt, and the 
entire structure is then separated from its periodic im- 
ages by 20 A of vacuum. In the isolated geometry, a 
dipole correction^ was applied in the center of the vac- 
uum region. 

The ferroelectric behavior of the PbTiC>3 films can be 
characterized by the set of rumpling parameters, Sz = 
•Zpb(Ti) — z o, which describes the cation-oxygen separa- 
tion along the polarization direction for each PbO and 
Ti02 layer in a given film. In continuous structures, 
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FIG. 1: Schematic of PbTi03 in the continuous (left) and 
isolated (right) geometries. Periodic slabs in the former are 
in contact to form a continuous lattice, while in the latter 
the slabs are separated by a large vacuum. The horizontal 
solid and dashed lines indicate the relaxed positions of the 
cations and oxygens, respectively, in the film surface normal 
direction. 



we expect that Sz will evolve as lpt increases, until the 
Pt thickness is sufficient to completely screen the two 
PbTiC>3 surfaces from each other. We expect the isolated 
geometry to behave similarly, with Sz changing until the 
electrodes are thick enough to prevent interactions be- 
tween each PbTiC>3 surface and the corresponding elec- 
trode/vacuum interfaces. Furthermore, we expect the 
converged Sz values to be equal in the two geometries. 

To determine the minimum electrode thickness for con- 
verged ferroelectric thin film properties, we consider the 
evolution of ferroelectric rumpling S function of 

electrode thickness Zpt. Table [I] gives the values for the 
continuous and isolated geometries, respectively, as a 
fraction of the bulk rumpling parameters. As the ta- 
ble shows, the ferroelectric behavior of the PbTiO,3 films 
with Pt electrodes is identical in the continuous and iso- 
lated geometries for structures with lp t > 7, with all 
atoms in the same positions to within <0.001 A. Simi- 
larly, the PbTiC>3 rumpling profile in the two geometries 
converge to a single structure with SrRu03 electrodes of 
nine or more layers. The table shows that PbTiC>3 films 
with Pt electrodes have polarization enhanced by w 15 % 
relative to bulk at the positive surface and decreased by 
slightly more at the negative surface, while the central 
layers of the film show bulk-like rumpling. The Sz values 
for films with SrRuC>3 electrodes indicate a suppression 
of the polarization in all layers, with rumpling parame- 
ters of w 20-35 % of the bulk values. 

In a grounded PbTiOs/metal system, electrons are 
able to flow through the electrode (s) from one ferroelec- 
tric/electrode interface to the other, regardless of ge- 
ometry. Consequently, an electric field cannot be sus- 



TABLE I: Rumpling parameters Sz for the continuous and 
isolated geometries as a fraction of the theoretical bulk rum- 
pling parameters, <5zbuik(TiO 2 )=0.51 and <5;zbuik(PbO)=0.87 
A. The first Ti0 2 layer corresponds to the bottom PbTiC>3 
surface in Fig. with the surface normal parallel to the 
polarization. The ferroelectric-paraelectric energy difference, 



AE 



is also shown, in units of eV. A negative 



value of AE indicates an energetically favorable ferroelectric 
state. 
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tained through the electrode, and therefore the electro- 
static potential must be flat through the electrode re- 
gion. Macroaveraging&^ii the electrostatic potential in 
the continuous geometry, as illustrated in Fig. shows 
that structures with electrodes of seven and nine atomic 
layers of Pt are indeed in short-circuit, with clear zero- 
slope regions through the center of the electrode. On 
the other hand, the structures with lp t < 7 do not ex- 
hibit rigorously equipotential regions. As the potential 
for Zpt = 3 illustrates, the Pt regions corresponding to 
each interface meet, leaving no equipotential region. In- 
creasing the Pt thickness to five layers allows separate Pt 
and interface regions to emerge, but a small electrostatic 
potential slope remains, demonstrating that Zpt = 5 is 
insufficient to completely prevent interactions between 
the periodic copies of the PbTiC>3 film. The small slope 
through the electrode for each Zp t is given in the left side 
of Table ^ quantitatively indicating the rapid decrease 
of the field with increasing electrode thickness. The zero 
slope through the electrode for films with lp t > 5 con- 
firms that these structures are in short-circuit. 
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FIG. 2: Layer-averaged electrostatic potentials, $(z), along (001) for structures in the a) continuous and b) isolated geometries. 
The curves for films with lp t = 3, 5, 7, and 9 are shown with black dotted, red dot-dashed, green dashed, and blue solid lines, 
respectively. The vertical black, red, green, and blue lines in a) bound a single repeat of the periodic supercell for lp t — 3, 5, 7, 
and 9, respectively; the lines on the left side of the ferroelectric film are positioned at the central Pt layer for each lpt (labeled 
in the corresponding color), and those on the right mark the same Pt layer in the next supercell. 



The electrostatic potentials of PbTi0 3 /Pt structures 
in the isolated geometry, plotted in Fig. |2t>, show a sim- 
ilar picture. As in the continuous geometry, it is not 
possible to clearly differentiate the electrode from the in- 
terfaces for Zp t = 3. In addition, five Pt layers are again 
not thick enough to eliminate fields through the electrode 
due to the polarization charge, as shown by the non-zero 
value of the slope given in Table [n] Most significantly, 
however, the potentials are again flat through the central 
region of the electrodes for the structures with Zp t > 5, 
suggesting that charge is successfully transferred between 
the two PbTi0 3 /Pt interfaces. The magnitude of the 
dipole correction, p corr , also reported in Table ITU further 
indicates that the system is in short-circuit; p C orr goes to 
zero for lp t > 5, as there is no net dipole through the 
Pt/PbTi0 3 /Pt structure. 

Fig. |21 also shows that the field through the ferroelec- 
tric slab, £fe, is very similar for all Zpt, regardless of 
the geometry. The magnitude of the field through the 
central unit cell (the middle three layers in Fig. QJ, re- 
ported in Table [HJ again indicates that the minimum 
electrode thickness for converged ferroelectric properties 
in both geometries is Zpt = 7. Interestingly, the evolu- 
tion of £fe with Pt thickness is dependent on the geom- 



etry, with EpE decreasing towards the converged value 
in the continuous geometry, and increasing to the con- 
verged value in the isolated geometry. The difference in 
the trend demonstrates the importance of choosing elec- 
trodes of sufficient thickness in studying the behavior of 
charge passivation in ferroelectric/metal capacitors. The 
identical converged values of -Efe, however, show that 
the two geometries are equivalent once this condition is 
satisfied. 

With SrRu0 3 electrodes, the field through the PbTi0 3 
is significantly larger than that found with Pt electrodes. 
Nevertheless, as Table |jj] shows, the evolution of the 
electrostatic potential as a function of electrode thick- 
ness is similar for both electrode materials. In the 
PbTi0 3 /SrRu0 3 films, the slope through the ferroelec- 
tric converges to a value of 0.053 eV/A and the field 
through the electrodes approaches zero as the SrRu0 3 
thickness is increased to ZsrRu0 3 =11, regardless of the 
film geometry. The larger minimum electrode thickness 
determined for the PbTi0 3 /SrRu0 3 films reflects the 
longer screening length of SrRu0 3 compared to Pt. 

The above data demonstrate that for both elec- 
trode materials, charge transfer occurs between the two 
PbTi0 3 /metal interfaces in the isolated geometry, de- 
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TABLE II: Magnitude of the electric fields J5 mota i and -Efe 
through the electrode and ferroelectric regions, respectively, 
as a function of electrode thickness for continuous and isolated 
geometries. The values for -E me tal in the isolated geometry 
show the average of the slope through the two electrodes. The 
magnitude of the dipole correction in the isolated geometry, 
p C orr, is also given, in units of eA. The electric fields, in units of 
V/A, are determined from the macroaveraged— electrostatic 
potentials by a linear fit to the central 4 A in the PbTiC>3 
films and the SrRuOa electrodes, and to the inner 1-2 layers 
in the Pt electrodes. The uncertainty in the electric fields is 
±0.0005 V/A. 
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spite the fact that the electrodes are not in direct con- 
tact. This is possible due to the nature of the computa- 
tion, which requires that the electrodes are provided with 
electrons from a single reservoir; i.e., there is a constant 
chemical potential of electrons. Consequently, the min- 
imum electrode thickness necessary to short-circuit the 
system is defined in the same manner as in the continuous 
geometry, as the thickness necessary to screen interac- 
tions between interfaces, in this case, the PbTiOa/metal 
and metal/vacuum interfaces sandwiching each electrode. 

Examination of the density of states (DOS) of the cen- 
tral metal layer offers further confirmation that short- 
circuit boundary conditions have been satisfied in both 
geometries. In grounded systems with electrodes of suf- 
ficient thickness to prevent interaction between periodic 
copies of the PbTiOa film, the electronic structure of the 
central metal layer should be identical to the correspond- 
ing layer in an (001) metal slab, as the effects of the ferro- 
electric film (or the vacuum) will be completely screened 
in this region. In Fig. El the d-band DOS of PbTi0 3 /Pt 
structures in the continuous and isolated geometries are 
compared to the DOS of the middle Pt layer in a 9-layer 
Pt (001) film with the same in-plane lattice constant. As 
the top panel of the figure shows, the d-band of the con- 
tinuous (solid curve) and isolated (dashed curve) struc- 
tures with Zp t — 3 differ significantly from the Pt (001) 
d-band (shaded). They also differ from each other, as ex- 
pected from the structural data shown in Table HJ As Zp t 
is increased, the continuous and isolated d-bands look 
increasingly similar to that of the Pt (001). Although 
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FIG. 3: d-band DOS of the central Pt layer in the electrodes 
of the continuous (solid red curve) and isolated (dashed blue 
curve) geometries. In both geometries, the Pt DOS in the 
center of the nine- layer electrodes is virtually identical to that 
in the center of a thick Pt (001) slab (shaded). 



there are still minor differences between the DOS of the 
7-layer structures and Pt (001), the isolated and continu- 
ous d-bands for the Ipt = 9 structures are both essentially 
indistinguishable from the Pt (001) d-band. 

The DOS of the central layer in the SrRuOs electrodes, 
not shown due to space constraints, also becomes very 
similar to the corresponding SrO or Ru02 layer in a thick 
SrRuOa slab as ZsrRu0 3 increases, displaying only small 
differences for ZsrRu0 3 > 9. 

In conclusion, we have shown that short-circuit bound- 
ary conditions are satisfied in ferroelectric PbTiOa /metal 
capacitor structures in two geometries when the electrode 
thickness is greater than five Pt layers or seven SrRuOa 
layers. We have shown that, for a given electrode ma- 
terial, the ferroelectric films in the continuous and iso- 
lated geometries converge to the same atomic structure, 
and that the potential in the ferroelectric and metal re- 
gions is identical for both geometries with electrodes of 
sufficient thickness, resulting in indistinguishable ferro- 
electric behavior. In such structures, the electrostatic 
potential has zero slope through the center of the metal 
electrode, demonstrating that electron transfer between 
the PbTiOa/metal interfaces is accomplished, via direct 
contact between the electrodes in the continuous geome- 
try or as a result of a single reservoir of electrons shared 
by the electrodes in the isolated geometry. An electronic 
structure analysis shows that the DOS of the central 
metal layer is almost identical to that in the correspond- 
ing metal slab, further demonstrating that these struc- 
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tures are in short-circuit. Our results demonstrate that 
ferroelectric capacitor structures can be accurately mod- 
eled in either geometry, allowing the possibility to utilize 
the benefits of both in further studies. 
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In fact, we find that perpendicular polarization is favored 
over in-plane polarization for PbTiOs in-plane lattice con- 
stants from 2.5% less than the equilibrium value up to 0.5% 
above the equilibrium lattice constant, at which point the 
in-plane polarization becomes more favorable. Since the 
Pt(100) lattice constant is only 0.38% larger than that of 
SrTiOa, we predict perpendicular polarization for PbTiC>3 
grown on Pt(100). The in-plane polarization has been stud- 
ied by I. Kornev et al. (Phys. Rev. Lett. 93, 196104 (2004) 
and by Meyer, Padilla, and Vanderbilt (Faraday Discus- 
sions 114, 395 (1999)). We here address the out-of-plane 
polarization while leaving the boundary conditions relating 
to in-plane polarization in a future study. 
We consider electrodes with an odd number of metal layers 
and a single middle layer, as required by the continuous 
geometry. 



